Decay on several sorts of heterogeneous 
centers: Special monodisperse approximation 
in the situation of strong unsymmetry. 4. 
Numerical results for the approximation of 
essential asymptotes 

V.Kurasov 
Victor.Kurasov@pobox.spbu.ru 

This manuscript directly continues 0], [Q, 0|. All definitions and formu- 
las have to be taken from The numerical results for comparison with the 
total monodisperse approximation has to be taken from p|. The numerical 
results for the floating monodisperse approximation can be taken from 0. 

1 Calculations 

We have to recall again the system of the condensation equations. It can be 
written in the following form 

G = [ exp{-G{x))6i{x){z - xfdx 
Jo 

9i = exp{—b / exp(— ^(a;))^^) 
Jo 

with a positive parameter b and have to estimate the error in 

roo 

N = / exp{-lG{x))dx 
Jo 

with some parameter /. 
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We shall solve this problem numerically and compare our result with the 
already formulated models. In the model of the total monodisperse approx- 
imation we get 

Na = / exp{—lG Aix))dx 
Jo 

where Ga is 

Ga = ^{1 - exp{-bD))x^ 
and the constant D is given by 

D = / exp(-x^/A)dx = 1.28 
Jo 

Numerical results are shown in 0. 

In the model of the floating monodisperse approximation we have to cal- 
culate the integral 

Nb = / exp{—lG Bix))dx 
Jo 

where Gb is 

Gb = t(1 - exp(-6 / exp(-x^/4:)dx))z^ 
b Jo 

Gb ^ t(1 - exp(-b(Q(D - z/A)z/A + Biz/A - D)D)))z^ 
b 

Numerical results are shown in [Q. 

It is very attractive to spread the approximation for the last integral at 
small z for all z (as it was done in the intermediate situation in Q when 
we solved the algebraic equation on the parameters of the spectrum (in the 
intermediate situation it is absolutely justified). Then we came to the third 
approximation 

Nc = / exp(~lGc(x))dx 
Jo 

where Gc is 

Gc ^ t(1 -exp(-te/4))z=^ 
b 

This approximation will be called as "approximation of essential asymp- 
totes". The real advantage of this approximation is the absence of the ex- 
ponential nonlinearity. When this approximation will be introduced into 
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equation on the parameters of the condensation process there will be no 
numerical difficulties to solve it. 

We have tried all mentioned approximations for b from 0.2 up to 5.2 with 
the step 0.2 and for / from 0.2 up to 5.2 with a step 0.2. We calculate the 
relative error in A^. The results are drawn in fig.l for Nc where the relative 
errors are marked by r^. 

We see that the relative errors of Nb and Nc are very small and practically 
the same. One can not find the difference between fig.l in |^ and fig.l here. 

The maximum of errors in Nb and Nc lies near / = 0. So, we have to 
analyse the situation with small values of /. It was done in fig. 2 for Nc- We 
see that we can not find the maximum error because it increases at small b. 
Then we have to calculate the situation with 6 = 0. Here we have to solve 
the following equation 

/•oo 

G= exp{-G{x)){z - x^dx 
Jo 

and to compare 

poo 

N = / exp{-lG)dx 
Jo 

POO 

Na = / exp{-lDz^)dz 
Jo 

POO 

Nb= / exp{-l{Q{z/4~ D)Dz^ + e{D - z/4)zy4))dz 
Jo 

roo 

Nc= / exp{-lzy4)dz 



with 





We can not put here / = directly. 

The results are shown in fig.3. One can see one curve with two wings. 
The upper wing corresponds to the error of A^i and the lower corresponds 
to the relative error in Nb and Nc- At / = these wings come together. 
We see that our hypothesis (the worst situation for the floating monodis- 
perse approximation takes place when the ffist type heterogeneous centers 
are unexhausted) is really true. 

The worst situation is when b is near zero and / lies also near zero. Here 
we can use the total monodisperse approximation to estimate the error. It 
is clear that the relative error in Na is greater than in Nb (not in Nc)- So, 
we can calculate ri, see that when b goes to zero it decreases (it is clear also 
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from the physical reasons) and estimate the error r2 at 6 = 0, / = by i?i 
calculated at small b and / = 0. Then one can see that it is small. 

An interesting problem is to see whether Nb and Nc are different or 
no. Earlier we can not see the difference. In fig. 4 one can see the ratio 
1^2/1^3 plotted at 6 = and can note that only for / ^ 0.01 ^ 0.02 one can 
see the small difference. It means that to see the difference between these 
approximations the ratio between the scale of the first type centers nucleation 
and the scale of the second type centers nucleation must be giant. Even at 
giant values the difference is small. 
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Fig.l 

The relative error of Nc drawn as the function of I and b. Parameter I 
goes from 0.2 up to 5.2 with a step 0.2. Parameter b goes from 0.2 up to 5.2 
with a step 0.2. 

One can see the maximum at small I and moderate b. One can not sepa- 
rate Nb and Nq according to fig.l in and fig.l. 
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Fig.2 



The relative error of Nq drawn as the function of I and b. Parameter I 
goes from 0.01 up to 0.11 with a step 0.01. Parameter b goes from 0.2 up to 
5.2 with a step 0.2. 

One can see the maximum at small I and small b. One can note that 
now the values of b corresponding to maximum of the relative errors become 
small. One can not separate Nb and Nq according to fig.2 in and fig.2 
here. 
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Fig.3 



The relative errors of Na, Nb and Nc drawn as the function of I atb — 0. 
Parameter I goes from 0.01 up to 5.01. 

One can see two wings which come together for b near 0. The upper wing 
corresponds to the relative error of Na- The lower wing corresponds to the 
relative errors of Nb and Nc- One can not separate them. 

One can see that ri decreases when b goes to and can estimate r2 by ri 
at small I. 
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Fig.4 



The ratio r2/r^ drawn as the function of I at b — 0. Parameter I goes 
from 0.01 up to 5.01. 

One can see the difference between r2 and only at very small values of 
b (only the first two points corresponding to I = 0.01 and I — 0.02^. So, N3 
can be also considered as a suitable approximation. 
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